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Abstract 

We refine a method for finding a canonical form of symmetry operators of arbitrary order for 
the Schrodinger eigenvalue equation = (A2 + V)^f = on any 2D Riemannian manifold, 
real or complex that admits a separation of variables in some orthogonal coordinate system. The 
flat space equations with potentials V = a(x + iy) k ~ 1 /{x — iy) k+1 in Cartesian coordinates, and 
V = ar 2 + /3/r 2 cos 2 k6 + 7/r 2 sin 2 kO (the Tremblay, Turbiner, and Winternitz system) in polar 
coordinates, have each been shown to be classically superintegrable for all rational numbers k. 
We apply the canonical operator method to give a constructive proof that each of these systems 
is also quantum superintegrable. We also develop the classical analog of the quantum canonical 
form for a symmetry. It is clear that our methods will generalize to other Hamiltonian systems. 



1 Introduction 

We consider an n-dimensional classical superintegrable system is an integrable Hamiltonian system 
that not only possesses n mutually commuting integrals, but in addition, the Hamiltonian Poisson- 
commutes with 2n — 1 functions on the phase space that are globally denned and polynomial in the 
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momenta. This notion can be extended to define a quantum superintegrable system to be a quantum 
Hamiltonian which is one of a set of n independent mutually commuting differential operators that 
commutes with a set of In — 1 independent differential operators of finite order. We restrict to classical 
systems of the form H = Y^7j=i 9^PiPj + ^ anc ^ corresponding quantum systems H = A n + V where 
V is related to V but in general is not equal to it, [H|2|. These systems, including the classical Kepler 
problem and the quantum hydrogen atom have great historical importance, due to their remarkable 
properties. They are exactly analytically solvable and in multiple ways. Superintegrable systems 
of 1st order, i.e., classical systems where the defining symmetries are first order in the momenta 
and quantum systems where the symmetries are first order partial differential operators, are directly 
related to Lie transformation groups and well understood. Superintegrable systems of 2nd order 
have been well studied and there is now a structure and classification theory [21 [3] H] [5]. The 
connection between 2nd order symmetries and separation of variables has been of crucial importance 
in finding examples and carrying out the classification [SI [71 [8]. However for 3rd and higher order 
superintegrable systems much less is known. In particular there have been relatively few examples 
and there is almost no structure theory, i.e., an understanding of the structure of the Poisson algebra 
generated by the classical symmetries or the algebra generated by the quantum symmetries and their 
commutators, and no classification theory. 

This situation has changed recently with the discovery of many more examples of classical (especially) 
and quantum superintegrable systems of order higher than two, j9j [TUl HH [12], [131 EE]- Also, the tool 
of coupling constant metamorphosis (Stackel transform) has been developed to map superintegrable 
systems of higher order on one Riemannian space to superintegrable systems of the same order and 
structure on a different Riemannian space [T51 dSl EL71 EIB1 EE] • In this paper we will concentrate on the 
case n = 2 where the number of independent symmetries is 3, including the classical Hamiltonian or 
quantum Schrodinger operator, respectively. In almost all of the families of new examples the second 
symmetry is of 2nd order and defines an orthogonal separable coordinate system for the classical 
Hamilton- Jacobi equation or the quantum Schrodinger equation. Only one defining symmetry is 
of higher order. We are particularly considering the classical example of Tremblay, Turbiner, and 
Winternitz [TTl [T2] where 

V = ar 2 + - + 1 

r 2 cos 2 k9 r 2 sin 2 k6 

in polar coordinates. Due to the separation in polar coordinates there are two commuting 2nd order 
constants of the motion. For certain rational values of the parameter k these authors found an addi- 
tional constant of the motion (usually of higher order), so that the system was superintegrable both 
in the classical and quantum sense. They conjectured and provided impressive evidence that these 
systems were classically and quantum superintegrable for all rational k. In [20] it was shown that, in 
fact all of the classical TTW systems were superintegrable. Quesne [21] used a structure developed 
by Dunkl to show that the TTW systems for k an odd integer were quantum superintegrable. As 
a bi-product of the tools developed in this paper we will give a constructive proof that the TTW 
system is quantum superintegrable for all rational k. However, our main contribution is a tool for 
the verification of classical and quantum superintegrability of higher order that can be applied to a 
variety of Hamiltonian systems. 

In Section 1 we review a construction of a canonical form for symmetry operators of all orders of a 
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time- independent Schrodinger equation that admits an orthogonal separation of variables 
Then in Section 2 we apply this tool to the flat space potential V = a(x + iy) k ~ 1 /(x — iy) k+l in 
Cartesian coordinates (separable in polar coordinates), that has recently been shown to be classically 
superintegrable for all rational k [20]. We demonstrate that it is also quantum superintegrable for 
all rational k. 

In Section 3 we give the analogous construction of a canonical form for symmetry operators of all or- 
ders of for a classical Hamiltonian system. We again treat the example V = a(x + iy) k ~ x j [x — iy) k+1 
and give a new proof that it is classically superintegrable for all rational k. The special interest here 
is the relation between the classical and quantum construction. 

In Section 4 we apply canonical operator method to the TTW case to demonstrate that it is quantum 
superintegrable for all rational k. In Section 5 we discuss our overall strategy and the prospects for 
exploitation and generalization of our methods. 



2 The canonical form for a symmetry operator 

We consider a Schrodinger equation on a 2D real or complex Riemannian manifold with Laplace- 
Beltrami operator A 2 and potential V: 

HV = (A 2 + V)V = E,V (1) 

that also admits an orthogonal separation of variables. If {ui,u 2 } is the orthogonal separable coor- 
dinate system the corresponding Schrodinger operator has the form 

H = L x = A 2 + V(u x , u 2 ) = 1 (d 2 Ui + d 2 U2 + Vl ( Ul ) + v 2 {u 2 )) . 

Jl\ U l) + J2{U 2 ) 

and, due to the separability, there is the second-order symmetry operator 



[d ui + viim)) - (d U2 + v 2 {u 2 )) 



/1N + /2W v Ul v " fi(ui) + f 2 (u 2 ) 
i.e., [L 2 ,H] = 0, and the operator identities 

h( Ul )H + L 2 = d 2 Ui + Vl ( Ul ), f 2 (u 2 )H -L 2 = d 2 U2 + v 2 (u 2 ). (2) 

We look for a partial differential operator L(H, L 2 , u\, u 2 ) that satisfies 

[H,L]=0. (3) 
We require that the symmetry operator take the standard form 

L = {A^( Ul ,u 2 )d ulU2 + Bi> k (u u u 2 )d Ul + C>' k ( Ul ,u 2 )d U2 + Di'\ Ul ,u 2 )) H j L k . (4) 



Note that if the formal operators (T4J) contained partial derivatives in U\ and u 2 of orders > 2 we could 
use the identities (j2j), recursively, and rearrange terms to achieve the unique standard form (J3J). 

Using operator identities and writing d u , = dj we have 

f 7/ f 11'- 

[d 1 ,H] = --^—H+—^—, [d 2 ,H] = --^-H + 



/1 + /2 /1 + /2' /1 + /2 /1 + /2 



/1 + /2 /1 + /2' 1 ' J /1 + /2 /1 + /2' 

Jl + J 2 Jl + J2 Jl + J2 

[H, F( Ul , u 2 )\ = —^—(F U1U1 + F U2U2 + 2F Ul d x + 2F U2 d 2 ). 

Jl + J2 

From these results and (jSJ) we obtain 

(/1 + A5 12 ] = 2A Ul (f 1 d 2 H + d 2 L 2 - Vl d 2 ) + 2A U2 (f 2 d 1 H - d 1 L 2 - v 2 d x ) 

+A(f 2 d 1 H + f[d 2 H - vfa - v[d 2 ) + (A U1U1 + A U2U2 )d l2 , 
(fi + f 2 )[H, Bdi] = B(f[H - v[) + 2B ux (hH + L 2 - Vl ) + (B UlUl + B U2U2 )d x + 2B U2 d 12} 
(fi + f 2 )[H, Cd 2 ] = C{f 2 H - v' 2 ) + 2C U2 (f 2 H -L 2 - v 2 ) + (C UlUl + C U2U2 )d 2 + 2C Ul d u , 
(fi + f 2 )[H, D] = D U1U1 + D U2U2 + 2£> ui di + 2D U2 d 2 . 

Thus we have 

(A(«l) + / 2 (w2))[#,^(Wl,«2)<9l2 + B(u l ,U 2 )d 1 + C(l4i,U 2 )3j + Z>(«i,tt a )] = 

(Atiwi + A U2U2 + 2B U2 + 2C Ul )d 12 + (B UlUl + B U2U2 — 2A U2 v 2 + 2D Ul — Av' 2 )di 
+ {2A U2 f 2 + Af 2 )d x H - 2A U2 d 1 L 2 + (C UlUl + C U2U2 - 2A Ul v x + 2D U2 - Av[)d 2 

+(2A u Ji + Af[)d 2 H + 2A Ul d 2 L 2 
+(D UlUl + D U2U2 - 2B Ul v x - 2C U2 v 2 - Bv[ - Cv' 2 ) 
+{2B u Jt + 2C U J 2 + Bf[ + Cf 2 )H + (2B U1 - 2C U2 )L 2 . 

The symmetry condition ([3]) is equivalent to the system of equations 

d n A j ' k + d 22 A^ k + 23 2 B j ' k + 2d 1 C 1 ' k = 0, (5) 

d n B j ' k + d 22 B j ' k - 2d 2 A j ' k v 2 + 2d x D^ k - A j < k v' 2 + 

(2d 2 A j ~ 1 ' k f 2 + A j - l > k f 2 ) - 2d 2 A 1 ' k ' 1 = 0, (6) 
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d u C j ' k + d 22 C^ k - 2d l A hk v 1 + 2d 2 D j > k - A j > k v[+ 

{2d 1 A j - 1 ' k h + A j ~ 1,k f[) + 2d 1 A 1 ' k ~ 1 = 0, (7) 

d u D j > k + d 22 D ] ' k - 2d 1 B j ' k v l - 23 2 C hk v 2 - B j ' k v[ - C^ k v' 2 (8) 
+(2d 1 B j ~ 1 > k f 1 + 2d 2 C j - 1 ' k f 2 + B j - l > k f[ + C j - l > k f 2 ) + (2d 1 B j ' k ~ 1 - 2d 2 C j ' k ~ 1 ) = 0. 

Note that condition (jlj) makes sense, at least formally, for infinite order differential equations. Indeed, 
one can consider H, L 2 as parameters in these equations. Then once L is expanded as a power 
series in these parameters, the terms are reordered so that the powers of the parameters are on the 
right, before they are replaced by explicit differential operators. Alternatively one can consider the 
operator L as acting on a simultaneous eigenbasis of the commuting operators H and L 2 , in which 
case the parameters are the eigenvalues. Of course (jH) is defined rigorously for finite order symmetry 
operators. 

In this view we can write 

L(H,L 2 ,ux,u 2 ) = A(ux,u 2 ,H,L 2 )d 12 + B(u 1 ,u 2 ,H,L 2 )d 1 

+ C(ui,u 2 ,H,L 2 )d 2 + D( Ul ,u 2 ,H,L 2 ), (9) 

and consider L as an at most second-order order differential operator in u±, u 2 that is analytic in the 
parameters H, L 2 . Then the above system of equations can be written in the more compact form 

A UlUl + A U2U2 + 2B U2 + 2C U1 = 0, (10) 

B U1U1 + B U2U2 - 2A U2 v 2 + 2D U1 - Av 2 + (2A U J 2 + Af' 2 )H - 2A U2 L 2 = 0, (11) 

C uxUl + C U2U2 - 2A UlVl + 2D U2 - Av[ + {2A U J 1 + Af[)H + 2A Ul L 2 = 0, (12) 

D U1U1 + D U2U2 - 2B Ul v x - 2C U2 v 2 - Bv[ - Cv' 2 (13) 

+{2B U J X + 2C U2 f 2 + Bf[ + Cf 2 )H + (2B U1 - 2C U2 )L 2 = 0. 

We can view ffTUj) as an equation for A,B,C and f|TT]) . f|T2l) as the defining equations for D Ul ,D U2 . 
Then L is L with the terms in H and L 2 interpeted as (jlj) and considered as partial differential 
operators. 

We can simplify this system by noting that there are two functions F(ui, u 2 , H, L 2 ), G(ui, u 2 , H, L 2 ) 
such that ffTUjl is satisfied by 

A = F, B — ~d 2 F - dtG, C = ~d x F + d 2 G, 

Then the integrability condition for ( TTTj) . ( TT2"|) is (with the shorthand djF = Fj, djiF = Fji, etc., for 
/•' and <7). 

2G 1222 + ^F 2222 + 2F 22 (v 2 - f 2 H + L 2 ) + 3F 2 (v' 2 - f' 2 H) + F(v% - ftH) = 
-2G 1112 + l -F lul + 2F n (v 1 - frH - L 2 ) + 3F X K - f{H) + F{v'{ - f»H), (14) 

5 



and equation f|T3|) becomes 

l -F lll2 + 2F l2 {v l -f 1 H) + F l {v' 2 -f 2 H)+ l -G llll + (15) 

2G 11 (v 1 -f 1 H-L 2 ) + G 1 (v[-f[H) = 

~F l222 - 2F 12 (v 2 - f 2 H) - F 2 (v[ - f[H) + X -G 2222 + 
2G 22 (v 2 - f 2 H + L 2 ) + G 2 (v' 2 - f' 2 H). 

We remark that any solution of (TH|) . (fl5|) with A, 5, C not identically corresponds to a symmetry 
operator that does not commute with L 2 , hence is algebraically independent of the symmetries H, L 2 . 



3 The potential V = a ^-L+i 

We consider the flat space Schrodinger operator 

H = d xx + d yy + V, V = a^±^, (16) 

where x, y are Cartesian coordinates. We have shown that the corresponding classical systems are 
superintegrable for all rational k, 



In the special case k = 3 we have explicitly established quantum superintegrability. Indeed, we 
obtained the symmetry operators 



Ki = (d x - idyf + _ . [-(zy + 3x)9 x + (3zy + x)d y ], 
K 2 = (xd y - yd x )(d x - id y f + — — — — [i(2y 2 - 3ixy - 3x 2 )d 2 x - (3iy + x)(iy + 3x)9 1 .<9 ?/ + 
i(2x 2 + 3ixy - 3y 2 )<9 2 - 2i(3iy + x)d x - 2(iy + 3x)d y - 8z] + ia 2 ^ + ^ 



(x - iyY 



f _ 7/ 2 I Q r 2A 

K 3 = (xd y - yd x f + 2iay± 1 " 



(x — iy) 3 



>2 -)2 ^,(a; + ^) ( 



H = dt + d, + a- 



v {x 2 + y 2 ) A 
with the commutation relations 

-2 



[#i, K 2 \ = 3iK(, [K lt K 3 ] = 6iK 2 - 9K lt 
[K 2 , K 3 ] = 3i{Kx, K 2 } + z(27 + 6a)^i + 9K 2 , 
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and the analogue of the constraint 

1 9 63 

-{K u K u K 3 } - 3K\ - i-{K u K 2 } + ( — + 3a)Kf - 3aH 3 = 0. 

All of these quantum systems separate in polar coordinates: 

U\ = R, u 2 = 6 ', x = e R cos9, y = e R sin9. 
The corresponding symmetry operator is 

—L 2 = K 2 = o e + ae 

Furthermore, 

H = e' 2R {dl - L 2 ) , 

and 

fx = e 2R , h = 0, v x = 0, v 2 = ae 2ike . 

We assume k = p/q for positive relatively prime integers p, q. Based on the known expressions for 
the classical higher order constants of the motion, derived in [], we look for an operator constant of 
the motion L, ^j, where 

F = J2^A^,H,L 2 )e^ R+ibke \ G = J2l3aA^H,L 2 )e 2 ^ R+ake \ (17) 

a, b a, b 

We require that there are only a finite number of nonzero terms in the sums and that the sums are 
of the form a = a + m, b = bo + n where m, n run over a subset of the non-negative integers. (Thus 
C ao ,b w iH be an analog of a lowest weight vector. Substituting all these expressions into equations 
piji ; ([1311 and equating coefficients of terms e 2 ( aR + lbk9 ) ; we obtain the matrix recursion 

,2 u2^2\ ( iakb a 2 + k 2 b 2 - L 2 \ ( A a ,b 



(2a-l)H 



-ikb 1 - a \ / Ai-i,6 
a i I B _i,6 



za — 1) \ / A a 



or, solving for 



(2a-l)g / a(L 2 + 3k 2 b 2 - a 2 ) 4z(l - a)aA;6 , , , 



J(a.b) V -*kb(L 2 + k 2 b 2 ) -(a-l)(L 2 + a 2 + k 2 b 
ak(2b - 1) / -£;o(L 2 - Po 2 + 3a 2 ) 4z(6 - l)ak 2 b 



J(a,b) V m(L 2 + a 2 ) fc(6-l)(L 2 + a 2 + fc 2 6 2 ) 



C a ,b-1 — 0, 



where 



J(a, b) = 2(a 2 - k 2 b 2 )((a - kbf - L 2 )((a + kbf - L 2 ). 



Consider first the case where p, q are both odd. We see from fTTHT) that we can choose the 2-tuple 
C-p/2, q /2 arbitrarily, i.e., it is not a consequence of a recursion. Thus we set a = —p/2, 6 = —q/2, 
so that <2g — k 2 b\ = 0. Further we set C a ^ = unless it can be computed explicitly from C ao ,&o by a 
sequence of recursions (ITU|) . 

Think of the elements C a ^ as laid out on an infinite grid, with rows labeled by a and columns by b. 
The value of C aQ+m ^ 0+n for m, ^ can be obtained via (IT3|) as the sum of the contributions from 
all regular paths that lead from gridpoint (a ,b ) to gridpoint a + m, 6 + n. A regular path is a 
connected sequence of vertical moves upward: (a — 1, b) — > (a, b) and horizontal moves to the right: 
(a, b — l)—> (a, b), in arbitrary order. All contributions from gridpoints below or to the left of (do, bo) 
are assumed zero. The contribution of a path to the value of C ao+m ^ 0+n is the ordered product of the 
contributions of the vertical and horizontal segments that make up the regular path. 

We have assumed first that p and q are odd, so that p/2, q/2 are half-integers. Now notice that 
any regular path connecting (a ,b ) = (—p/2, —q/2) to some gridpoint (1/2,6) necessarily contains 
a vertical segment (—1/2,6) — > (1/2, 6) and contributes the factor 0, so that the contribution of the 
path to the sum is 0. Thus, necessarily, C\nj> = 0. Similarly, any regular path connecting (a , 6q) = 
(—p/2,— q/2) to some gridpoint (a, 1/2) necessarily contains an horizontal segment (6,-1/2) — > 
(6, 1/2) and contributes the factor 0, so that the contribution of the path to the sum is 0. Thus, 
C a ,i/2 = 0. We conclude that the only possible nonzero 2-tuples are those in the grid (—p/2 + 
m, —q/2 + n) where < m < p/2, < n < q/2, and these terms are uniquely determined by the 
choice of C- p /2,- q /2- We get polynomial constants of the motion by taking the terms in C- p /2,- g /2 to 
be suitable finite products of the form 

]J[J(a,b)}, 

a,b 

to cancel the denominator terms in the the expressions for C ao + m ,b +n which come from recursion (TL9]h 
Thus we have constructed a 2-parameter family of finite order constants of the motion. It is a simple 
exercise to show that the commutators of these symmetries with L 2 are nonzero, so that the system 
is operator superintegrable. (Note that this last fact follows also from our original construction of 
the symmetries L. We must have A = B = C = unless L is functionally independent of H and 
L 2 .) 

If k = —p/q with p, q both odd, the same construction works with (a , 6 ) = (—p/2, —q/2). 

The case k = 2 s p/q with s > 1 and p, q relatively prime odd integers requires a modified analysis. 
Now we set ao = — 2 s ~ 1 p, 6 = —q/2, so that dp — /c 2 6q = 0, a Q is an integer and, as before, bo is half 
integer. Further we set C a ^ = unless it can be computed explicitly from C ao ,b by a sequence of 
recursions (TH?]) . As before, the value of C ao+m ^ Q+n for m, ^ can be obtained via (IT9|) as the sum of 
the contributions from all regular paths that lead from gridpoint (do, bo) to gridpoint ao + m, bo + n. 
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Now notice that any vertical segment connecting a gridpoint (—1, 6) to gridpoint (0, 6) maps C_ x ~ b to 




i.e., to a 2-vector with upper component 0. Further, if this segment is followed by the horizontal 
segment connecting (0,6) to (0,6 + 1) the upper component of the 2-vector will remain 0. Thus all 
regular paths that lead from gridpoint (ao, bo) to any gridpoint (0, b) on row a = will produce a 
2-vector of the form 

Next, note that any vertical segment (0,6) — > (1,6) will map a special 2-vector fl20|) to the zero 
vector. This means that C a ^ = for all integers a > 1. Just as before, C a ^ = for all half-integers 
6 > 1/2. Thus the only possible nonzero 2-tuples are those in the grid (— 2 s ~ 1 p + m, —q/2 + n) where 
< m < 2 s ~ 1 p, < n < q/2, and these terms are uniquely determined by the choice of C_ p /2- q /2- 
We get polynomial constants of the motion by taking the terms in C- p / 2 - q /2 to be suitable finite 
products of the form 

Y[[J(a,b)], 

a,b 

to cancel the denominator terms in the the expressions for C ao _|_ m! fc + n . Thus we have again con- 
structed a 2-parameter family of finite order constants of the motion and the system is operator 
superintegrable. 

It is easy to extend these arguments to the cases k = —2 s p/q and k = ±p/2 s q where p, q are relatively 
prime odd integers. Thus the system (TT6l) is operator superintegrable for all rational k. 



4 The classical analog 

Here we first describe the classical analog of our infinite order symmetry operator construction and 
then apply it to the same example as in the previous section. We construct constants of the motion 
of all orders for the Hamiltonian system 

2 

-H = 9 3 %Pk + V = E (21) 
j,k=l 

that admits a separation of variables. If {«i,W2} defines an orthogonal additive separable coordinate 
system for the Hamilton- Jacobi equation in some Riemannian space, the corresponding Hamiltonian 
system has the form 
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and, due to the separability, there is the second-order constant of the motion 

r /2O2) / 2 , / x\ hM ( 2 , / n\ 

£a = AW + AM (P '" + l( 0) " /,(«,) + « (? " + 1!< " 2)) ' 

i.e., {£2,%} = 0, where {•, •} is the usual Poisson bracket, and we have phase space identities 

h{ Ul )H + £ 2 = p 2 Ul + «i(u0, / 2 (u 2 )ft - £ 2 = < + t^(« 2 ). (22) 

We look for a constant of the motion £(H, £2, ui, u 2 ), i.e., a function on the phase space that satisfies 

{U,£} = 0. (23) 
We require that the constant of the motion take the standard form 

£ = (A j ' k (ui,u 2 )p Ul p U2 + B^( Ul , u 2 ) Pui + C^( Ul , u 2 ) Pu2 + Df> k ( Ul , u 2 )) W£\. (24) 

Note that if the formal symmetries ()24p contained polynomial terms in p Ul or p U2 of orders > 2 we 
could use the identities f )22|) . recursively, and rearrange terms to achieve the unique standard form 



We find that the symmetry condition (1231) is equivalent to the system of equations 

d 2 B j ' k + <9iC J ' fe = 0, (25) 

- 2d 2 A j ' k v 2 + 2d l D j ' k - A j ' k v' 2 + 2d 2 A j ~ 1 ' k f 2 + A j ~ l > k f 2 - 2d 2 A j > k - 1 = 0, (26) 

- 2«9iA i '% + 2d 2 D j ' k - A j ' k v[ + 2d l A j ~ 1 ' k f 1 + A j ~ l ' k f[ + 2d 1 A j > k ~ 1 = 0, (27) 

- 2d 1 B j > k v 1 - 2d 2 C j ' k v 2 - B hk v[ - C hk v' 2 + 2d 1 B j ~ 1 > k f 1 (28) 
+2d 2 C J - 1 ' k f 2 + B'-^fi + C j - l ' k & + 2d l B j ' k ~ 1 - 2d 2 C j ' k - 1 = 0. 

Note that condition (|24|) makes sense, at least formally, for infinite order constants of the motion, 
and one can consider H, £ 2 as parameters in these equations. 

In this view we can write 

£CH,£ 2 ,u 1 ,u 2 ) = A(u 1 ,u 2 ,n,£ 2 )pip 2 + B(u 1 ,u 2 ,'H,£ 2 )p 1 

+ C( Ul ,u 2 ,U,£ 2 )p 2 + D(u u u 2 ,H,£ 2 ), (29) 

and consider £ as an at most second-order constant of the motion that is analytic in the parameters 
'H,£2- Then the above system of equations can be written in the more compact form 

B U2 + C U1 = 0, (30) 

- 2A U2 v 2 + 2D U1 - Av' 2 + {2A U2 f 2 + Af 2 )H - 2A U2 L 2 = 0, (31) 
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- 2A ulVl + 2D m - Av[ + {2A U J X + Af[)H + 2A Ul L 2 = 0, (32) 
- 2B UlVl - 2C U2 v 2 - Bv[ - Cv' 2 (33) 
+(2S U1 /! + 2C U J 2 + Bf[ + Cf 2 )H + (2B U1 - 2C U2 )L 2 = 0. 
We can view (!30|) as an equation for B, C and (I3~TT) . (|32|) as the defining equations for D U1 , D U2 . 

We can simplify this system, and easily compare it to the operator system, by noting that there are 
two functions F(ui, u 2 , H, C 2 ), G(ui, u 2 , %, C 2 ) such that (1301) is satisfied by 

A — F, B = -dxG, C = d 2 G, 

Then the integrability condition for fl3TT) . fl32|) is (with the shorthand djF = Fj, djgF = Fji, etc., for 
FandG), 

2F 22 {v 2 - f 2 U + C 2 ) + 3F 2 {v' 2 - f 2 H) + F(v' 2 ' - ftU) = 

2F 11 {v 1 - hU - C 2 ) + 3F 1 (v[ - f[H) + F(v'l - f{%), (34) 
and equation fl33|) becomes 

2G n (v 1 -f 1 n-JC 2 ) + G 1 (v' l -f[n) = (35) 
2G 22 (v 2 - j 2 U + C 2 ) + G 2 (v' 2 - f' 2 U). 



Now we use this classical construction to study the flat space Hamiltonian system 

(x + iy) k ~ 1 
(x - iy) 



H=pl+pl + V, V = a )*Z (36) 



where x, y are Cartesian coordinates. We have already shown that this system is superintegrable for 
all rational k, |2T)1. 



All of these classical systems separate in polar coordinates: 

Ui — R, u 2 — 9, x = e R cos9, y = e R sin9, 
with corresponding constants of the motion 

-C 2 =p 2 + ae 2M . 

Furthermore, 

n = e~ 2R ( P % - c 2 ) , 

and 

fi = e 2R , f 2 = 0, v x = 0, v 2 = ae 2lke . 

We assume k = p/q for relatively prime integers p, q. Based on the known expressions for the classical 
higher order constants of the motion, derived in [20J, we look for a standard form constant of the 
motion £, (|2"91) . where 

F = Y,A a , b (a,n,C2)e 2{aR+lbke \ G = £ Bafi (a, % C 2 )e 2 ^ bkd \ (37) 

a,b a,b 
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We require that there are only a finite number of nonzero terms in the sums and that the sums are of 
the form a = a Q +m, b = bo+n where m, n run over a subset of the non-negative integers. Substituting 
all these expressions into equations (1551) and equating coefficients of terms e 2 ( aR + lbk9 ) ) we obtain 
the matrix recursion 

2C 2{ »- - *V) ( £ ) + (2a -1)E [I a °_ 1 ) ( ) (38) 

Although this system of equations is much simpler than the corresponding operator equations (jl8p . 
ffl9|) . it shares essential features with them so that the details of the proof of superintegrability are 
essentially unchanged. As before we set 




Consider first the case where p, q are both odd and positive. We see from (138]) that we can choose 
the 2-tuple C_ p / 2 ,<j/2 arbitrarily. Thus we set a = —p/2, b = —q/2, so that Oq — k 2 b\ = and we set 
C a ,b = unless it can be computed explicitly from C ao ,£>o by a sequence of recursions (13"5]) . 

The value of C ao+mt b + n for m, ^ can be obtained via (|38|) as the sum of the contributions from 
all regular paths that lead from gridpoint (a ,b ) to gridpoint a + m,b + n. Since p/2, q/2 are 
half-integers, any regular path connecting (a ,b ) = (—p/2, —q/2) to some gridpoint (1/2,6) nec- 
essarily contains a vertical segment (—1/2,6) — > (1/2, b) and contributes the factor 0, so that the 
contribution of the path to the sum is 0. Thus, necessarily, Ci/2,6 = 0. Similarly, any regular path 
connecting (ao,&o) — (—p/2, —q/2) to some gridpoint (a, 1/2) necessarily contains an horizontal seg- 
ment (b, —1/2) — > (b, 1/2) and contributes the factor 0, so that the contribution of the path to the 
sum is 0. Thus, C a i/2 = 0. We conclude that the only possible nonzero 2-tuples are those in the grid 
(— p/2 + m, —q/2 + n) where < m < p/2, < n < q/2, and these terms are uniquely determined by 
the choice of C_p/ 2> _ g / 2 . We get polynomial constants of the motion by taking the terms in C_ p / 2j - g /2 
to be suitable powers of £ 2 to cancel the denominator terms in the the expressions for C ao+m> b 0+n 
which come from recursion 038 p . Thus we have constructed a 2-parameter family of finite order 
constants of the motion. It is easy to show that the Poisson brackets of these symmetries with £ 2 
are nonzero, so that the system is classically superintegrable. There is a special simplification here 
in that the recursion fl38|) decouples into separate equations for A a ,b and for B a ,b- 

If k = —p/q with p, q both odd, the same construction works with (a , b ) = (—p/2, —q/2). 

For case k = 2 s p/q with s > 1 and p, q relatively prime odd positive integers, we set ao = — 2 s ~ 1 p, 
b Q = —q/2, so that — k 2 b\ = 0, a is an integer and, as before, b is half integer. Again we set 
C a ,b = unless it can be computed explicitly from C ao ,fco by a sequence of recursions ( 13"S|) . As before, 
the value of C ao + m ,b +n for m, ^ can be obtained via (138]) as the sum of the contributions from 
all regular paths that lead from gridpoint (ao, bo) to gridpoint a + m, b + n. However, any vertical 
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segment connecting a gridpoint (—1,6) to gridpoint (0, 6) maps C_ x g to 



r - - , 





0,6 



i.e., to a 2-vector with upper component 0. If this segment is followed by the horizontal segment 
connecting (0, 6) to (0, 6 + 1) the upper component of the 2-vector will remain 0, so all regular paths 
that lead from gridpoint (do, &o) to any gridpoint (0, 6) on row a = will produce a 2-vector of the 
form 

Note that any vertical segment (0,6) — > (1,6) will map a special 2-vector fl39|) to the zero vector. 
Thus C a> b = for all integers a > 1, and as before, C 0) & = for all half-integers 6 > 1/2. Thus the 
only possible nonzero 2-tuples are those in the grid (— 2 s ~~ l p + m, —q/2 + n) where < m < 2 s ~ 1 p, 
< n < q/2, and these terms are uniquely determined by the choice of C_ p / 2j _g/ 2 . We get polynomial 
constants of the motion by taking the terms in C_ p / 2 ^ q j 2 to be suitable powers of C 2 to cancel the 
denominator terms in the the expressions for C ao +m,b +n- We have again constructed a 2-parameter 
family of finite order constants of the motion and the system is classically superintegrable. There is 
a significant simplification here, due to the decoupling of ( |38|) into separate equations for A a ,b and 
for B a £- If we choose B a0t b = then all vectors Co,& — 0, so that row can be removed from the grid. 

Again it is easy to extend these arguments to the cases k = —2 s p/q and k = ±p/2 s q where p,q 
are relatively prime odd integers. Thus we have a new proof that the system (|36|) is classically 
superintegrable for all rational k, and have clarified the relation between the classical and operator 
symmetries for this system. 

Example: For system ( |36|) with k = 3 the recurrence relations can be solved easily to give 
F = C 2 e- 3R - 3id + ^He- R -™ } G = C 2 e- m ~ 3W + ^e^ 3 * 6 ', 

so 

A = F, B = 3£ 2 e- 3 ^ 3<e +^e^ R ^, C = -MC 2 e~ 3R ~ 3W ~He~ R ~ 3W , D = l -C 2 (4C 2 +3e 2R H)e- 3R - 3W . 

This gives us the canonical forms for the 3rd and 4th order constants of the motion /C 1; /C 2 as given 
in 



Indeed, Ap R p e + D = f/Ci, Bp R + Cp e = \K, 2 



5 The TTW system 

Now we apply our constructions to the TTW system, [TTJ H2] . Here, 

u x = R,u 2 = 6, fx = e 2R , f 2 = 0,vx = ae m , (40) 
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_ 7 . 2( 7 + g) 2( 7 - g) cos(2£;fl) 

U2 cos 2 (^) sin 2 (^) sin 2 (2^) sin 2 (2A;fl) 

Based on the results of [20J for the classical case, we postulate expansions of F, G in finite series 

F = A aj b, c E a ,b,c(R, 0), G = B a>bjC E aj b tC (R, 9), (41) 

a,b,c a,b,c 

E a ,b,o = e 2aR sm b (2k9), E aM = e 2aR sm b (2k6) cos(2k9). 

The sum is taken over terms of the form a = a$ + m, b = bo + n, and c = 0, 1, where m, n are integers. 
The point (ao, bo) could in principle be any point in M 2 , however, for reasons discussed below, we will 
take a to be a positive integer and b to be a negative integer. 

Taking coefficients with respect to the basis ()4ip in each of equation ffl~4j) and (jl5|) gives recurrence 
relations for these coefficients. For example, the coefficient of e 2aii sin b-2 (2/c#) cos(2fc(9) in equation 
( Jlij) gives the equation 

8bk 2 {b - 1) (L 2 - 2{k 2 {b 2 + 1) + 7 + 0)) A aM + 32ak 3 b (b 2 - l) B a , b+lfi 

+ 8k 2 (b 2 - 1) (6 2 A; 2 + 26A; 2 + 2(3 + 2 7 ) A a , b+2il + 16A; 2 (b 2 - l) (7 - /3) A a , b+2i0 

- 8A; 2 (2b -l)(b- 1) (7 - [3) A aA o + 8 (a 2 - A; 2 (6 - l) 2 ) (L 2 - a 2 - fc 2 (6 - l) 2 ) A a ^ x 

+ AHa (2a - 1) A a _ li6 _ 2jl - 8aa (a - 1) A a _ 2j6 _ 2il + 32aA; (6 - 1) (a 2 - k 2 (b - l) 2 ) S aj6 _ 1)0 = 0. 

The shifts in the indices of A and B are integers and so we can view this as an equation on a 
two-dimensional lattice with integer spacings. While the shifts in the indices are of integer size, we 
haven't required that the indices themselves be integers, although they may be integers in particular 
examples. 

Taking the coefficient of e 2aR sin fe (2fc#) in equation (j!4p and the coefficients of e 2aR sm b ~ 1 (2k9) cos(2k9) 
and e 2aR sin 6 ^ 1 (2k9) in equation f[T5"j) give a further three recurrence relations. At a general point 
in the lattice there are 4 coefficients, and these 4 equations will be shown to be independent. The 
equations are linear and homogeneous and so there must be some points where the independence of 
the equations breaks down and allows at least one coefficient to be arbitrarily chosen. 

The different powers of sin(2k8) used in obtaining these equations have been chosen as a matter of 
convenience after many experiments conducted using the computer algebra package Maple. 

All four recurrence relations are of a similar complexity, but rather than write them out separately, 
we will combine them into a matrix recurrence relation by defining 

/ A a,b,0 \ 

n B a ,b-1,0 
<^a,b — a 

Ai,6-2,1 
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Figure 1: Points contributing to the recurrence relation are marked with large dot (•). The large 
dot in the upper left corner corresponds to the position (a, b). 

We can now write the 4 recurrence relations in matrix form as 

= M a ^C a ,b + M a ^-2C a ,b~2 + M a ^-iC a ,b-i + M a ^~6C a ,b-6 

+ Ma-i^Ca-Xj) + M a _i ) 6 +2 Co-l,6+2 + M a _2,bC a -2,b + M a ^ 2 ,b+2C a-2,6+2, (42) 

where each Mjj is a 4 x 4 matrix given below. It is useful to visualize the the set of points in the 
lattice which enter into this recurrence for a given choice of (a, b). These are represented in figured] 
in which the upper left corner is the point (a, b). From this it is clear that when M a ^ is nonsingular, 
the value of C ai & can be uniquely determined from the 8 points to its right and below. In that case, 

Ca,b = ~M a j ) (^M a ,b-2Ca,b-2 + Mi,b-4Ca,b-4 + M afe _ 6 C ai fe_ 6 

+ M a -i t bCa-l,b + M a ^i^+ 2 Ca-l,b+2 + M a _ 2 ,bC a -2,b + ^a-2,6+2 Ca-2,fe+2 ) • (43) 



This allows us to construct an iterative procedure that calculates the values of C a ^ at points in the 
lattice using only other points where the values of Cjj are already known. Since the point (a, b) 
corresponds to the top left corner of the collection of points in Figure [TJ this process will calculate 
the coefficients in a sequence that moves from right to left and bottom to top. Note the matrices 
corresponding to the right hand 'corners' of the set of points in Figure [1] (M aj b +6 and M a _ 2 ,b+2) are 
singular and so could not be used in the same way, while the matrix corresponding to the bottom 
left corner is generically not singular, however it has properties that we will use for another purpose. 



M, 



a,b 



( 8A; 2 (26-i)(6-l)(/3- 7 ) 32ak (b - 1) (a 2 - k 2 (b - i f) 

8 (a 2 - k 2 b 2 ) (La -a 2 - b 2 k 2 ) 

8abk (a 2 - k 2 b 2 ) 

V 8ak (2b - 1) 08 - 7) -8 (a 2 - k 2 (b - l) 2 ) (L 2 - a 2 - k 2 (b - l) 2 ) 

(a 2 -k 2 (b-l) 2 )(L 2 -a 2 -k 2 (b-l) 2 ) \ 

-32abk (a 2 - k 2 b 2 ) 

-8 (a 2 - A; 2 6 2 ) (L 2 -a 2 - b 2 k 2 ) 

-8ak(b-l)(a 2 -k 2 (b-l) 2 ) 8bk 2 (2b - 1) (f3 - 7) ) 
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/ 16k 2 (6 2 - 1) (7 - p) 32ak 3 b {b 2 - 1) 

8 {b + 1) k 2 ((6 + 2)(L 2 - 2k 2 {b 2 + 2b + 2)) - 2b((3 + 7)) 

8ak (6 + 1) (k 2 (b + 2) + 2(f3 + 7)) 8k 2 (b + 1)(26 + 1)( 7 - P) 

\ Wak (b + 1) (7 - p) 8bk 2 (b+l)(2(k 2 (b 2 + l) + /3 + 1 )-L 2 ) 



8bk 2 {b-l){L 2 -2{k 2 (b 2 + l)+(3 + 1 )) \ 

8bk 2 (2b + 1) (7 - (3) 32ak {b + 1) (a 2 - 2k 2 (b + 2) 2 ) 

8aA; (26 + 1) (/3 - 7) 8 (6 + 1) A; 2 (2(6 + 2)(/3 + 7) + 26A; 2 (6 2 + 2b + 2) - 6L 2 ) 

8a6A;(a 2 -2(A; 2 (6 2 + l) + /3 + 7)) 8A; 2 {b + 1) (46 + 3) (7 - p) ) 



M, 



a,fe+4 








/ 

8A; 2 (b + 3) (b + 1) (A; 2 (6 + 4) (6 + 2) + 2(f3 + 7)) 

16fc 2 (6 + 3)(6 + l)(/3-7) 

-8A; 2 (6 + 3) (6 + 1) (6A; 2 (6 + 2) + 2(/3 + 7)) 







8k 2 (6 - 1) (6 + 1) (6A; 2 (6 + 2) + 2(/5 + 7)) 

8A; 2 (46 + 5) (6 + 1) (/3 - 7) 32aA; 3 (6 + 3) (6 + 2) (6 + 1) 

16ak (6 + 1) (7 - P) -8k 2 (6 + 3) (6 + 1) (bk 2 (b + 2) + 2(/3 + 7)) 
8ak (6 + 1) (6A; 2 (6 + 2) + 2(f3 + 7)) lQk 2 (6 + 3) (6 + 1) (f3 - 7) 



/ 



/ 



Ma. 



,6+6 







\ 



16k 2 (6 + 3) (6+1) (7-/3) 





V 










/ 



/ 

4#a (2a - 1) 
-4/J6A; (2a - 1) 
\ 









AHa (2a - 1) 





-AH (2a - 1) (a - 1) 4if A; (6 - 1) (2a - 1) 



\ 


-4/7 (2a - 1) (a - 1) 




,6+2 



/ \ 





\ -AHbk (2a -1)0/ 
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M, 



a-2,b 















— 8aa (a - 


- 1) 







— 8aa (0 — 


1) 


o 






o 






o 


8akb (a — 


1) 
















8a (a — 1) (a 





8a (a - 


-1) 


(a — 


2) 


-8a/c(6- 1) 


Co- 





















oN 








































































U 





8akb (a — 1) 









/ 



We are interested in finding a solution to the recurrence relation that gives F and G as finite sums 
and hence we seek solutions that are confined to a finite rectangle in the lattice. Since the equations 
for the Aij^ arid -Bjj^ are linear and homogeneous, they always admit the trivial solution and so we 
need to demonstrate that a nonzero solution can be found. Our approach is as follows. 

For a finite solution, there must be a lowest nonzero row and in that row, a rightmost nonzero 
element. Label this rightmost point in the bottom row as (do, bo). Since all elements to the right and 
below this point are zero, M a0t b must be singular, otherwise we could use (143]) to show that C ao ,b 
must vanish, contradicting its definition. 

Since 

det(M 0j6 ) = -4096 (a 2 - k 2 b 2 ) 2 (a 2 - k 2 (b - l) 2 ) 2 (L 2 - (a + k(b - l)) 2 ) (L 2 - (a - k(b - l)) 2 ) x 

x (L 2 - (a + kbf) (L 2 - (a - kbf) (44) 

we must choose our starting point so that either a 2 = k 2 b 2 or a 2 = k 2 (b — l) 2 , that is, if k = p/q 
with p and q a pair of relatively prime positive integers, we can choose (a , b ) to be one of (r]p, rjq), 
(— rjp, rjq), (r]p,r]q + 1) or (rjp, —rjq + 1) for any real number rj. At all of these points, the rank of 
M aQt b is 2 and hence at these points we can choose 4 — Rank(M aj fc) = 2 components of C a ,b to be 
arbitrary parameters. 

In order to have a finite solution, we must eventually reach a point in the lattice beyond which all 
entries vanish or can be chosen to vanish. Examining the matrices defining the recurrence, we see 
that it may be possible to achieve this a left hand boundary due to the many terms with factors such 
as 6 — 1, 6+1, 6 + 3 and on an upper boundary because of factors of a and a — 1. For this reason, 
we will now take a to be — p and 6 to be q or q + 1 and examine how the cut offs on the left and 
top occur. 

As the recurrence relations fj4*2|) and ( 143|) only involve shifts of multiples of two units in second index, 
it is easy to verify that all entries in columns that are an odd number of steps away from column 6o 
must vanish or can be chosen to vanish. Furthermore, we have two candidates for 6 , q and g + 1. We 
will choose 6q to be which ever of these is odd. We can then assume that even numbered columns 
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have only vanishing entries and can traverse the lattice in steps of two to the left starting from 
column 6 . 

We now work our way across row ao starting from column 60 taking steps of two units to the left. 
At the first point, (a,b) = (a ,6 ), M a>b has rank 2 and so the components of C a ^ depend linearly 
on two arbitrary parameters. At other points in the bottom row, M a ^ is nonsingular (unless we 
reach b = — 6 ) and so we can solve for C aQ ^. At the points when b = 1,-1,-3, this takes a 
special form. Note that in the bottom row, all lower points have vanishing Cij and so we need only 
consider contributions from the points (a, b + 2), (a, 6 + 4) and a, 6 + 6). We will initially assume that 
q ^ 1,2, 3, 4, 5, 6 so that M_ P) i, M_ Pi _ 1 , M_ Pj _ 3 and M_ Pi _ 5 are all nonsingular. This is not essential, 
however the argument is simpler in this case. 

First consider 6 = 1. The form of the matrices giving contributions from points to the right are the 
first two matrices in Table [TJ It is clear from these that the third component of C a0) i must vanish. 
Next consider 6 = —1. The only nonzero matrix elements occur in column 3 of M ( ^ ) 1 _ 1 M a0jl and so 

Coo,-! = M-^M^C^ + M-;_ 1 M ao , 3 C ao ,3 + M-^M^C^s = 0. 
A similar calculation shows that C a0i _3 = C a0i _5 = and hence C ao j = for all j ' < — 1. 

Next we repeat the process for the row above, that is, row a + 1 starting from the right hand end, 
and then again for row ao + 2 and so. The argument showing that all Cjj vanish for j < —1 is 
essentially the same as for row ao except there are a few extra terms to consider since the elements 
in the row below are no longer all zero. 

To see how the cut off occurs at the top, start at the right hand end of the row, that is in position 
(0,6o). All elements to the right are zero and so the only contributions to Cq^ come from below, 
that is from (— l,6o) and (— 2,6o). It is clear from the corresponding matrices in Table [2] that the 
first and third components of Co,b are zero. Stepping across the row in step of two to the left, it is 
easy to check that this is maintained for all elements of this row. 

Next consider row 1 . From the form of the matrices given for a = 1 in Table [2] it is clear that the only 
the first and third components of Coj can contributed to any Ciy. However, since these components 
have already been shown to vanish we conclude that Cij = for all j. 

The last step is to consider row 2. As for row 1, when a = 2 the form of the matrices given for a = 2 
in Table |2] clearly shows that only the first and third components of Cqj can contributed to any C^j'- 
As these components have already been shown to vanish we conclude that C2J = for all j. 

Since we have two completely zero rows, it is now clear that Cij = for all z > 1. 

The above argument needs modification to see that the left hand cut off can be achieved when 
treating the bottom row for q = 1,2, 3, 4, 5 or 6 as M a ^ will be singular in one of the columns 6 = — 1, 
—3 or —5. However, it is a simple matter to use the original matrix recurrence relations (j4"2j) to check 
that the same conclusions can be reached in each of these cases, that is, the third component of C a(h \ 
vanishes and each entry of C ao j for j < 1 is either required to vanish or can be chosen to vanish. 

Example 1: It is well known that the TTW system is quantum superintegrable in the case k = 2, 



18 



1: 



/ * * * * \ 

* * * * 


* * * * 



M'lM a , b+2 



( * \ 

* 



\ * ) 

M~lM a , b+6 



( * * \ 

* * * 

0*0 

\ * * / 



( \ 

0*0 



\ / 

M-lM a ^ b+2 



6= -1: 



/ * 

* 

* 

\ * 

M-b M a,b + 2 



0\ 






/ \ 





\ / 



/ \ 

* 

0*0 

\ / 

M" b 1 M a _ li6+2 
M- 1 M a _ 2)6+2 



b= -3: 



/ * \ 

0*0 

0*0 

\ * ) 



( \ 





\ / 



b=-5: 



( * 

* 

* 

\ * 



\ 



0/ 



Table 1: Matrices giving contributions near the left boundary. A '*' represents a nonzero entry. 
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a = 0: 



/ * * \ 

* * 

* * 

\0 * * / 

M~lM aM . 2 



( * \ 



0*0 

\ / 



( \ 

* * * * 



\ * * J 



/ \ 

* 



\ J 

M~lM a . 1M2 
M-}M a _ 2 , h+2 



a = 1: 



/ * \ 

* * 

* * 
\ * / 



( \ 

0*0 
0*0 

\ o o o o ) 

M~}M a ^ h+2 



/ \ 





\ o o o o / 

m:Im, 



0,6 
L a,b 



a-2,b 



MzlM a . 2 u 2 



2: 



/ * \ 

* * 

* * 
\ * J 

M-lM a ^ b 



( \ 

0*0 

0*0 

\ ) 

M- 1 M a _ 2i6+2 



Table 2: Matrices giving contributions near the top boundary. A '*' represents a nonzero entry. 
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(p — 2, q = 1), [TT], [20] . The generating operators, expressed in Cartesian coordinates are 

{x 2 + y 2 ) 



H = d 2 + 3 2 + a(x^ + y z 



(xd y -yd x ) 2 + 4(3 



(3- 



(x 2 + y 2 ) 



x 2 y 2 



(x 2 — y 2 ) 2 



x z — y 
+ 7 



2\2 



+ 7" 



4x 2 y 2 



> 4 + y 4 ) ,0,7 
4x V P 2 ' 



L = (d 2 — d 2 ) 2 + (2ax 2 + 2/3 



(x 2 + y 2 ) (x 2 — y 2 



(x 2 — y 2 ) 2 



-7" 



2x 2 y 2 



-iaxy + 



7 



l(3xy 



(x 2 - y 2 ) 



r)^ + (2m/ 2 + 2/3 



(x 2 + y 2 ) (x 2 — y 2 ) 
2x 2 y 2 



7 



+ (2ax - ^)d x + (2m/ - + a 2 (x 2 - y 2 ) 2 + 



(x 2 — y 2 ) 2 
P 2 



+ 7- 



,2\ 2 



(x 2 — y 2 ) 2 



Y{x 2 -y 2 Y 



16x 4 y 4 



8a/3 



x 2 y 2 



[x 2 - y 2 ) 2 2x 2 y 2 71 2x 4 2y* 1 



By expressing the 4th order symmetry L in polar coordinates and converting to canonical form we 
can read off the functions A, B, C, D and then determine F, G and the nonzero expansion coefficients. 
The results are 



C. 



2,1 



/ _44 - 4L 2 \ 

4( 7 - /3) 


\ -28 - 8L 2 J 



, C-i, 



( -2H\ 





Co A 



( o \ 




\2a ) 



It is easy to check that these terms satisfy all our recurrence relations. 



Example 2: The nonzero vectors in the solution to the recurrence for the k = 1/3 quantum TTW 
system are C-\ t \, C-1,3, Co,3 and Co,i- The solution to the recurrence depends linearly on two 
parameters that can be taken to be y4._2,i,o and B -2,0,1- To obtain solutions for A, B, C and D that 
are polynomial in H and L 2 , we must choose the free parameters so as to cancel any demoninators. 
With the choice 



B_ 



2,0,1 



36#(18L 2 + 13) and A_ 2 ,\,o = 8(81L 2 2 + 765L 2 + 274), 



we find that L is a 6th order symmetry and the expressions for A, B, C and D are given below 
written in terms of the U\ = R and u 2 = 6 coordinates. 



A 



5 (81L 2 2 + 765L 2 + 274) e~ 2ui sin 3 Qu 2 J - 144( 7 - 0) (9L 2 + 26) e~ 2ui sin QuJ cos Qu 2 
+ 6 (6 (23 + 9L 2 ) (7 + /3) - 81( 7 - /3) 2 - 81L 2 2 - 765L 2 - 274) sin Q« 2 J e~ 2ui 
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B = +40 (8 + 81L 2 2 + 135L 2 ) e~ 2ui sin 2 Qw 2 J cos Qu 2 J 
+ 48 (101 + 144L 2 ) ( 7 - /3)e" 2ui sin 2 (|u 2 



(12 (315L 2 + 229) (7 + 0) - 2754( 7 - 0) 2 - 80 - 810L 2 2 - 1350L 2 ) e~ 2ui cos Qu 2 J 
+ 12 (162(7 - /3)( 7 + /3) - (450L 2 + 319) (7 - 0)) e~ 2ui 

C = -40 (8 + 8LL 2 2 + 135L 2 ) e~ 2ui sin 3 \%i 2 \ + 144 (27L 2 + 8) (7 - 0) cos Qit 2 J sin Q« 2 J e" 2 " 1 
+ 6 (6 (-27L 2 - 5) (7 + /3) + 81( 7 - /3) 2 + 405L 2 2 + 40 + 675L 2 ) e~ 2ui sin Q« 2 J 
- 72# (182* + 13) sin 3 Q« 2 J + 129622( 7 - /9) sin Q« 2 J cos Qu 2 
+ 5422 (-6(7 + /3) + 18L 2 + 13) sin Qu 2 

D = -(l2(81L 2 2 + 423L 2 + 40)(7 + /3)-162(9L 2 + 8)( 7 -/3) 2 )cosQM 2 ')e- 2ui 

+ 12 (9L 2 + 40) (9L 2 + 1) (7 - 0) cos (^u 2 j e~ 2ui + 2 (81L 2 2 + 765L 2 + 274) L 2 cos (2u 2 ) e~ 2ui 
(36 (-27L 2 - 23) (7 - /3)( 7 + 0) + 162( 7 - /3) 3 + 6 (81L 2 2 + 441L 2 + 40) (7 - 0)) e~ 2ui 
+ 81 (-2 (3L 2 + 7) (7 + 0) + 9(7 - 0) 2 ) if cos Qw 2 J + 486 (L 2 + 2) 22( 7 - /3) cos ^~u 2 J 
+ (81L 2 2 + 441L 2 + 40) 22 cos (2u 2 ) 

These give the operator L ([9]). In order to construct the symmetry operator in standard form (J4J), 
the 22 and L 2 , which have been treated as parameters throughout the calculation, must be moved 
to the right. For example, after expanding C, the coefficient of HL 2 is 

1296 sin 3 ( |w 2 ] + 972 sin ( \u 2 



,3 ) V3 

and so this contributes the term 

•1296 sin 3 ( %u 2 ] +972 sin ( ~u 2 ] ] -^-HL 



3 / \3 J J du 2 



2 



to the differential operator L, in which 22 and L 2 are now treated as differential operators. We have 
used Maple to verify that the this operator does in fact commute with the k = 1/3 quantum TTW 
Hamiltonian. 
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6 Discussion 



Key to our method for proof of superintegrability is the canonical form for symmetry operators 
of all orders. It enables us to replace the computation of the commutator of H with operators of 
arbitrary high order by verification of equations ( 114)) and ( II 5p . In these equations H and L 2 can be 
treated as parameters until the very last step when the canonical form is reinterpreted as an operator. 
Since ()14p and (fl~5l) are linear and homogeneous in F and G the solutions of these equations form 
a vector space. There are, of course, many solutions but most are not polynomials in H,L 2 . To 
prove superintegrability we have to find a nontrivial solution F(ui, u 2 , H, L 2 ), G(ui, u 2 , H, L 2 ) that 
has polynomial dependence on H, L 2 . If there is one such solution, there will be an infinite number of 
others, since any polynomial function of a finite symmetry is a finite symmetry, as is the commutator 
of L 2 with a finite symmetry. To prove superintegrability we need find only one such solution. We 
choose the simplest ansatz that leads to success. The method we employ will not necessarily lead to 
the symmetry operator of lowest order. 

Our strategy is to postulate a set of basis functions and to expand F and G in terms of them. 
The basis has to be chosen so that (TT4")) . (Tl5j) reduce to a set of recurrence relations between the 
coefficients of the basis functions. We will succeed if we can find some nonzero solution of these 
recurrences such that only a finite number of the coefficients are nonzero. The coefficients will then 
be rational functions of H, L 2 , but arbitrary up to a scale factor K(H, L). We choose K such that all 
coefficients become polynomials in H, L 2 , and then we are done! We have used two different methods 
to solve the recurrences in our two examples. The template method for the more complicated TTW 
problem is more general and its step-by-step evaluation of the expansion coefficients probably makes 
it the preferred tool to treat additional examples. 

How can one determine an appropriate set of basis functions? In the examples appearing in this paper 
we used the known expressions for the corresponding classical superintegrable systems as computed 
in [20] to determine the basis for the classical constants of the motion, and then used the same basis 
for the quantum system. This worked although the classical expansion coefficients differed from the 
quantum coefficients, as would be expected. It is clear that the methods of this paper will apply 
to many Hamiltonian systems, but each system will have its own peculiarities. Also, the canonical 
form for symmetry operators can clearly be extended to higher dimensions in the cases where the 
separable coordinates are of the subgroup type as treated in |14j . Of particular interest is the relation 
between the classical constants of the motion and the quantum symmetries. We intend to pursue 
these lines of inquiry. 
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